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FOREWORD 


This document was prepared by Northrop Services Inc., Huntsville, Alabama, 
for the Systems Dynamics Laboratory of the Science and Engineering Directorate 
of Marshall Space Flight Center, Huntsville, Alabama. Technical coordination 
was maintained through Mr. Mario H. Rheinfurth of the Dynamics and Trajectory 
Branch (ED15) , and Mr. Pat Vallely and Ms. Alberta King of the Servo Mechanisms 
and Systems Stability Branch (ED14) . The material presented herein is the 
result of work performed under Contract NAS8-29627 and is submitted as an 
interim report in response to the requirements of the said contract. 
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ABSTRACT 


This report documents both an analytical stability study and a digital 
simulation response study of two connected rigid bodies. Relative rotation of 
the bodies at the connection is allowed, thereby providing a model suitable 
for studying system stability and response during a soft-dock regime. Pro- 
visions are made for a docking poi;t axes alignment torque and a despin torque 
capability for encountering spinning payloads. Although the stability analysis 
is based on linearized equations, the digital simulation is based on nonlinear 
models . 
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Section I 

INTRODUCTION 


An important baseline mission of the Space Tug is the recovery of satel- 
lites at or below synchronous orbital altitudes for return to the Space Shuttle. 
Such recovery type missions for the Tug are to include the recovery of mal- 
functioning satellites. A probable malfunction is partial or complete loss of 
vehicle attitude control. The Tug must be able to dock and successfully cap- 
tu/fi with a high probability on the first attempt. 

In order to assess the feasibility of recovering a spinning satellite 
that is in a stable spin, but lacks the capability to despin, a study is 
made of the stability and response of the post-dock configuration. That is, 
the vehicles are initially assumed to already be in a soft-dock regime wherein 
capture or initial latch has been achieved but final latch or hard-dock has 
not. This system configuration is modeled as two rigid bodies connected at a 
single point about which relative rotation is permitted. The spin axis of the 
satellite (payload) is restricted to lie in the vicinity of the docking port 
axis. The behavior of this system model is then assessed relative to mass/ 
inertia characteristics, vehicle coupling through the docking mechanism (joint 
alignment torque) , despin torquer limits for despinning the payload, and ideal 
3-axis attitude control of the Tug. 

A dual approach is taken to the above problem. First, an analytical 
stability study is made based primarily on a linearized model thereby enabling 
the utilization of such linear stability techniques as root locus. This 
analysis makes up Part I of this report and is presented in Sections II 
through VI. In Section II, an introduction to the analytical stability study 
is given. Secondly, a simulation response study based on a nonlinear system 
model is presented in Part II (or Section VII) . Finally, in Section VIII 
a summary of results for both the analytical and simulation study are presented 
along with conclusions and some discussion. 
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Section II 

INTRODUCTION TO THE ANALYTICAL STABILITY ANALYSIS 


2.1 GENERAL 

Upon considering two rigid bodies connected at a single point and each 

spinning primarily along its body x-axis, it is only natural to be concerned 

with the relative motion between the x-axes of the two bodies . If an initial 

misalignment between these two axes continually grows with increasing time, 

then the coupled spinning body system is said to be unstable. In this context, 

"stability" is taken to mean the stability of the misalignment angle between 

the two-body spinning axes and it is readily seen that this misalignment angle 

is directly related to the y and z-directional rotations of each body, that 

is 0. , 6. , 0_ , and 0„ . In Figure 2-1, the Euler angles 0. , 0. , and 0. 

(i = lj2) are illustrated for a 3-2-1 (z-y-x) rotation sequence. It is seen 

that the position of the x^-axis is located relative to the inertial reference 

frame X, Y, Z by two angles 0. and 0. . 

“ xy xz 



Figure 2-1 . ROTATIONAL COORDINATES OF BODY i 
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As with most cases in nature, this is a nonlinear system with state vari- 
ables intermingled, thereby complicating the stability analysis. Although 
such theories as Popov's criterion and Liapunov's second method exist for non- 
linear system stability analysis, they are either not suitable for this specific 
problem or difficult to apply. Consequently, a more popular and easier method 
known as "linearization" is adopted even though it is often times less powerful 
than Liapunov's second method. 


The idea is to study the linear effect of perturbation around an equilib- 
rium operation point in the state space. If, in addition to the previously 
stated assumption that the two bodies are spinning primarily along their body 
x-axes , the y and z-directional rotations are assumed small, then the equilib- 
rium operation point becomes obvious, namely that for which each body is 
spinning only about its x-axis and the y and z-directional rotations are null 
(perfect alignment between the body fixed reference frames) . Thus 


ly 


= 0 


and 


Iz 


2y 


= ®2z = 


ly 


= 9 


Iz 


2y 


= 6 


2z 


= 0 . 


By dropping out all of the second and higher order terms involving 0. ,0- , 

• * • • 

02 y, ®2z’ ^ly’ ^Iz’ ^2y* ^2z* appropriate linearized system of equations 

is obtained (see Section III). A study of these linearized equations of motion 
would then hopefully reveal the stability of the original nonlinear system. 


The theoretical foundation of this linear analysis is Liapunov's first 
method (refs. 1, 2), which can be summarized as follows: 

• If the linearized system is asymptotically stable, the (asymptotic) 
stability of the nonlinear system is assured. 

• If the linearized system is unstable, the nonlinear system is also 
unstable. 

• If the linearized system is stable but not asymptotically stable, the 
stability of the nonlinear system will be determined by higher order 
terms . 


2-2 



NORTHROP SERVICES* INC. 


TR-1400 


For a free-spinning two-body system (one for which no external environ- 
mental torques or stabilization devices are considered) , it is quite possible 
that the linearized model cannot determine nonlinear stability. However, the 
objective of the study is to learn how to drive the misalignment angles close 
to zero in a reasonably short time, that is, to control the system to be 
asymptotically stable. Since a spring and damper are to be used to align and 
stabilize the x-axes of the spinning system, it is expected that the linearized 
model will prove to be both a powerful and convenient tool for analysis. 

There is still another question which remains to be answered, that con- 
cerning the "region of stability." The stability region is defined as the 
maximum region in the state space for which any initial condition will result 
in a trajectory that approaches the equilibrium point with increasing time. 

It should be noted that, in general, nonlinear asymptotic stability is only a 
local property rather than global, although linear stability is always global. 
Unfortunately, Liapunov’s first method cannot provide any information on the 
stability region and this question must then be investigated separately. 

2,2 VEHICLE ATTITUDE CONTROL, ALIGNMENT, AND DESPIN TORQUE OPTIONS 

For both demonstrative and practical purposes, a number of operational 
conditions relative to vehicle attitude control, alignment, and despin torque 
are to be considered in the analytical stability study. Relative to vehicle 
attitude control (A), there are three options which can be imposed on body 1. 
These are given in Table 2-1, For the A-2 control option, the desired conno- 
tation of "x-axis only" attitude control is implied by the condition 6^^ = 0. 
Although the A-3 control option indicates no attitude control, this does not 
preclude the possibility of body 1 having an initial spin rate identical to 
that of body 2. In addition, there are four options available relative to an 
alignment and/or despin torque (L) . These are presented in Table 2-2. Together, 
Tables 2-1 and 2-2 form a matrix of 12 operational conditions to consider. It 
should be noted that L-1 and L-2 are open- loop options whereas L-3 and L-4 are 
closed loop options . 
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Table 2-1. BODY 1 IDEAL ATTITUDE CONTROL OPTIONS 


A-1 

Body 1 has ideal "3-axis" attitude control 

A-2 

Body 1 has ideal "x-axis only" attitude control 

A-3 

Body 1 has "no" attitude control 


Table 2-2. ALIGNMENT AND/OR DESPIN TORQUE OPTIONS 


L-1 

Without alignment torque/without despin torque 

L-2 

Without alignment torque/with despin torque 

L-3 

With alignment torque/without despin torque 

L-4 

With alignment torque/with despin torque 


In the following analytical stability study, a simple but practical model 
for the despin torque is used. For the A-1 and A-2 control options, a constant 
x-directional despin torque (L^^) is assumed to be acting as long as the x-direc- 
tional relative spin rate between the two bodies (Aw) is not zero (see Figure 
2 - 2 ) . 



k 

W A 1 


^ LlOJ 

-1 



Figure 2-2. DESPIN TORQUE MODEL 
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2.3 NOMINAL PARAMETERS FOR THE TWO-BODY MODEL 

In the following study, analysis and design is based on a "nominal model" 
with variant models examined as necessary. Table 2-3 presents the parameters 
defining the configuration of the "nominal" two-body model. 


Table 2-3. NOMINAL TWO-BODY MODEL 



Body 1 

Body 2 


(Tug) 

(Spinning Satel li te) 

Shape 

”’t 


Mass 

m-j = 3222 Kg 

m 2 = 4500 Kg 

Inertia 

= 8418 Kg-m^ 

B^ = = 31000 Kg-m^ 

A 2 = 43400 Kg-m^ 
82 = 02 = 11732 Kg-m^ 

■ ■ 

Joint Vector 

d] = 5 m 

d, = d, =0 
ly Iz 

d2x “ -2-3 
d2y 


Each body is assumed to be symmetrical about the body x-axis with respect to 
inertia resulting in an inertia matrix of the form 

I. = 

X 

The bodies are assumed to be homogenous thereby requiring the center of mass 
of each body to lie on the x-axis. In addition, the joint connecting the two 
bodies is assumed to lie on the x-axis of each body resulting in a joint vecto 
(vector from the center of mass to the joint) for each body of the form 


A. 

3 . 


C. 0 

X 

0 c, 

X 


, i = 1, 2 
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, i = 1, 2 . 


The nominal value of the initial spin rate of body 2 is 0.1 rad/sec along the 
X2~axis . 
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Section III 

LINEARIZED EQUATIONS OF MOTION AND STATE VARIABLE CANONICAL FORM 


3.1 LINEARIZED EQUATIONS OF MOTION 

In Appendix A, a derivation is given for the nonlinear equations of 
motion for two rigid bodies connected at a single point. In Appendix B, 
these nonlinear equations of motion are simplified according to the following 
assumptions : 

• No external environmental forces or torques are acting on the bodies. 

• Partial linearization; The rotations about the y and z-axes of both 
bodies are assumed small as well as the products of these angles 
and their rates. 

• Both bodies are symmetrical about the body x-axis with respect 
to inertia. 

• The center of mass of each body is located on its respective body 
x-axis . 

• The joint or connection between the bodies is located on the x-axis 
of each body. 


Based on these assumptions, the resulting two-body equations of motion are 
given by 

A, •• A, /N 

J 6 + K 0 = L (3-1) 

with 


CJ 

0 

-K 

0 


C* 

0 

-K 


-K 

0 


0 

-K 

0 

c* 


K 


-Ai6, 

1 lx 
0 
0 


A.0, 

1 lx 


0 

0 


0 

0 


■^2®2x 


^2®2x 
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“ — 




— 

®iy 


^ly 


ly 




A 

i-i 

Iz 


. 9 = 

®2y 

, and L = 

A 

— T 



^2y 

®2z 


®2z 


-^2z 


where 


L. = L. cos 0. - sin 9. 

ly ly XX iz ix 


and 


L, = L. sin 6. + L. cos 0. 

iz xy XX iz XX 

C*, = C, + m*d? 

X X X 


i = 1, 2 


K = m*d^d2 


m 


* = 





e, ,0. ,0. 

ix xy xz 


0 < » 0 - » 0 . 
ix xy xz 

0 . , 0 . , 0 . 

XX xy xz 


- principal x-axis inertia of body 1, 2 

- principal y and z-axis inertia o£ body 1, 2 

- distance along the x-axis from the body 1 , 2 center 
of mass to the joint between the bodies 

- joint torque about the y and z axes of body i, 
i = 1, ,2 

- mass of body 1,2 

- body rotations about the x, y, and z axes of body 
i, i = 1, 2 

- angular velocities of body i, i = 1, 2 

- angular accelerations of body i, i = 1, 2 
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3,2 JOINT TORQUE/ALIGNMENT AND DESPIN 

The joint torque, L, as expressed in inertial space is given by 


L 




^ 2^-2 


where and are the coordinate transformation matrices resolving body 1 
and 2 coordinates, respectively, into the inertial reference frame (see 
Appendix B for the derivation of and T^) . To a first order approximation, 
the joint torque expressed in the inertial frame is then given by 



L. 

XX 


L, 

XX 

L. - 

— X 

L. 

xy 

= 

L. 0. + L, COS0. - L. sin0. 

XX xz xy XX xz xx 


L. 

xz 


-L. 0. + L. sin0 . + L, cos0. 

XX xy xy xx xz xx 


Comparison of this equation to those defined in subsection 3.1 for and 

L. results in 
iz 


L. = L - L e. 
ly y X iz 


and <3-2) 


L. - L + I 6. 

xz ,z X xy 


which will prove useful in subsection 3.3. 

The joint torque results from two separate sources, namely an alignment 
torque or spring-damper feedback that will tend to align the x-axes of the 
two bodies and a despin torque which will be applied to the x-axis of one of 
the bodies for despinning. In general then, the alignment torque (L ) and the 

A 

despin torque resolved in body 1 will have the form 



0 


^Dix 

^Al' 

.^Aiy 

V' 

0 


^Aiz 


0 
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i- = I^A 5=0 


Based on the despin torque model presented in subsection 2.2, the despin 
torque acting on body i can be writeen as 

-Di " ’ i = 1, 2 


where 

I - magnitude of the despin torque (see Figure 2-1) 

i. - X-axis unit vector of body i 
i 

then 


L^.= JiT. i. 

— Dx X X 


or 




T 


1 

ill 



T.„, 

= 1 

0. 

121 

xz 

T 


-0 

131 


ly 


, i = 1, 2 


(3-3) 


To damp the scissoring-type angular motion between the bodies, an align- 
ment or spring-damper feedback torque of the form 

L = k V + A V 
“A ~ ~ 

is provided where k and A are spring and damper coefficients and is an error 
vector based on the amount of angular misalignment between the x-axes of tlie 
two bodies. That is. 


V = Ij X 
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which, when resolved into the inertial reference frame becomes 



If two new coordinates are defined as 
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or 


-A 
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0 


— - 

0 

= k 

n 

+ X 

n 






4 


c 


J=A- 


L. 


0 

Ax 



% 

= 

kn + Xn 

^Az 

L j 


k? + xt 


(3-4) 


3.3 REDUCED EQUATIONS OF MOTION AND STATE VARIABLE CANONICAL FORM 

In Appendix C, the simplified two-body equations of motion as given by 
equation (3-1) are even further reduced by making a change of variables. These 
reduced two-body equations of motion are 




Y 


Z 

-f 

h 


C 





0 


-Pi 


0 -p. 


?2 0 


0 -?3 0 


-P 4 0 


<Hy - ^2^ 

+ ^2> 

-<"2 


(3-5) 


with 


Y = + (l-a)02y 

Z - 0.63^ + (l-a)02^ 

^ = V 


' = ®2z - 

Pi ' % ®lx + ^ H'> 


P 2 '■^2 h ®lx ■ *2 ®2^ 

P3 - (a^ A3 63^ - X3 A^ 02^)/A 
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P 4 - (4 h ®lx 4 ^2 ^ 


a = a^/(a^ + a^) 

a. = C* - , i = 1, 2 

X X 

and 

A = C* C* - . 

It should be emphasized that the variables n and ? indicate the relative mis- 
alignment between the body 1 and 2 y-axes and z-axes respectively. The vari- 
ables Y and Z are analogous to translational coordinates of the center of 
mass of the two body system. 

With the use of equations (3-2) and (3-3) , the forcing term of the equa- 
tions of motion (3-5) reduces to 


Ac/ (a^ + a^) 

- An/ (a^ + a^) 

-[(a^ + a^) “ ^1 ^ 

-[(a^ + a^) + a^ n A]/A 


(3-6) 


if the despin torque is applied in body 1 and 


- Ac/(a^ + a^) 

- An/ (a^ + a^) 


(3-7) 


-[(a^ + a^) L^y + C ^-l/A 
-[(a^ + + a^ n <?<]/A 

if the despin torque is applied in body 2. It is easily seen that when there 
is no despin torque, the forcing term reduces to 


w 
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0 

0 


^ 2 ^ % 
+ " 2 > 



With equations (3-6) and (3-7), equation (3-5) can now be written as 




g = ^'/(aj^ + 

and 

h = (a^ + a^)/h 
where the subscript 

i = 1, when the despin torque is applied in body 1 
i = 2, when the despin torque is applied in body 2 . 


Consequently, the state variable canonical form of the linearized model in 


terms of phase variable format is, in general, given by 
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This equation coupled with the following x-directional equations of motion (for 
rotation about the x-axes of the bodies) , 



J 
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Section IV 

A-1 CONTROL OPTION - BODY 1 HAS THREE AXIS ATTITUDE CONTROL 


4.1 REDUCED EQUATIONS OF MOTION 

For the A-1 control option, whereby body 1 has ideal three axis attitude 
control, the constraint equations 


lx 


= 0 


ly 


= =0 
Iz 


can be used to further simplify the two-body equations of motion. Applying 
these conditions to equation (3-1) results in a set of four equations 


-K 9 
-K 0 

c* e 
c* 0 


2y 

2z 

2y 

2z 



+ A, 



- ^2 «2x 


2z 


2y 




(4-1) 


the last two of which are sufficient for studying the behavior of body 2 
relative to body 1. With the use of equations (3-2) and (3-3), then 


= L. - J10„ + £0, 

2y Ay 2z Iz 




if the despin torque is applied in body 1 and 


Lo = La 
2y Ay 


^2z ^Az 


if the despin torque is applied in body 2. Since 0^^^ and are both constants, 
then there is no loss of generality if they are assumed equal to zero. There- 
fore, the last two equations of (4-1) can be rewritten as 
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0 -CO 0 


CO 0 0 


or in state variable canonical form 








0 0 1 0 0 


0 0 0 1 0 


0 0 -CO 0„ 

X 2y 


-£ , 0 CO 

X 


oJL“: 


-i/c* 


-i/c* 


where 


C* 2x " 


*■1 C * 


(when the despin torque is applied to body 1) 


£ 2=0 (when the despin torque is applied to body 2). 

• • 

With the above assumption that 6^^^ = 0^^ = 02y = ® 2 y " then the alignment 
torque given by equation (3-4) reduces to 

V " ^ ®2y ^ ®2y 

'*2z + ^ ®2z 

and the state variable canonical form can be still further reduced to 


1 0 
0 1 


-k 

£. 

X 

-A 

-CO 

A 

A 

A 

A 

-£. 

-k 

CJ 

-A 
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4.2 LINEAR MODEL OPEN-LOOP BEHAVIOR, DECOMPOSITION, AND INTERPRETATION 

In this subsection, the open loop behavior (that is, without the spring/ 
damper alignment torque) of the linear model will be investigated first without 
a despin torque (in subsection 4.2d) and secondly with a despin torque (in 
subsection 4.4.2). In both cases, body 1 is assumed to have ideal 3-axis 
attitude control. Numerical calculations will be based on the nominal model 
as presented in subsection 2.3. 


4.2.1 Open-Loop Behavior Without the Despin Torque (A-1, L-1) 

For the case of the linear model with neither an alignment nor a despin 
torque, k = X = £ = 0 and equation (4-2) reduces to 



(D = e„ (a constant) . 
o 2x0 

Computation shows that there are four open-loop poles on the imaginary axis : 



4-3 


NORTHROP SERVICES, INC. 


TR-1400 


= 0.2003 j 
P, = -0.2003 j • 

As was mentibned in subsection 2.1, no definite conclusion about the total 
system stability can be made by Liapunov's first method if the linear model 
cannot be shown to be unstable or asymptotically Stable. The appearance of 
the multiple zero poles, however, further motivates the study at hand. For if 
both of the zero poles appear in the or 62 ^ process, the system would be 
linearly unstable and hence, nonlinearly unstable. If they do not, then the 
linear analysis cannot reveal the overall stability. 


By rearranging the order of state variables in equation (4-2) , the state- 
space description becomes 


?2y 

.®2y 


®2z 





1 

0 

-b 


0 i 0 
b I 0 
0 I 0 
“H — 


s 

% 

®2z 


0 


2 z 


with the output variable 


•» ■ [‘ 


"] 


!2y 

!2y 

®2z 


0 


2 z 


From the partitioned matrix, it is readily seen that the 02 ^ process is 

unobservable (references 1, 3) from the 0 process and that the 0 „ process 

zy zy 

can be described by a state space of dimension three. That is, the subsystem 


..2y 

®2y 




0 1 0 
0 0 b 
0 -b 0 


!2y 

^2y 


ana 
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can be separated from the original system. With this reduced state-space 
description, which can be shown to be completely observable, it can be seen 
that there are three modes in the 0_ process, namely 



= bj 

P3 = -bj 


which implies only linear stability. Since similar results can be shown for 
the 02 ^ process, no immediate conclusions about overall stability can be made 
based on a linear analysis. 


Due to the special structure of the system, however, a study of the 0 


2y 


and © 2 ^ process is warranted 


Again by rearranging the order of state vari- 
ables, the state-space description becomes 


- 


1 — 


” • 



0 b 0 0 


K 

..2y 


1 


.2y 

So 


-b 0 I 0 0 


0„ 

2z 


1 


2z 

0 


1 0 1 0 0 



2y 




2y 

• 


1 

0 1 1 0 0 


0„ 

^ 2z. 




L 2zJ 


with the output 


< 

• CD 

L-. 


1 0 j 0 0 

2 y 




II 

0 1 1 0 0 


0 

0 

0 

0 


2y 

2z 

2y 

2z 



After separating this equation into the subsystems 
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and 





- 



r. 

-1 




0 

b 


0. 


2y 






2y 


”” 


-b 

0_ 




L 2z_ 



u 



L ZzJ 



r 

-1 



-I 




1 0 


®2y 


e„ 


0 1 


6 



— 2z_ 


L 

J 



2z 

• 



• • 

it is readily seen that the 02 ^ and processes are governed by two modes: 
bj and -b j , which means sustained oscillation in the precession rates. 

4.2.2 Open-Loop Behavior With the Despin Torque 

The affect on the stability of a single rigid body upon application of a 
despin torque has been discussed in subsection 3. 3. 2. 2 of reference 4. The 
case for the above linear model without an alignment torque but with a constant 
despin torque is discussed in the following. 


Considering first, the x-directional equation of motion, equation (3-10) 
^2 ® 2 x "" "^ 2 x " ® 2 x ° 


then 


and 


e = - ^ 

2 x 


2 x 


-T— t + ui , 0 <t<tT, 

o — — D 


where the despin time, t^^, is given by 




2 o 


and 03^ is the value of 02 ^ at t = 0. The appearance of 62 ^ in the coefficients 
of the y and Z directional equations of motion (equation 4-2) makes the system 
time-varying. Although the trajectory of the system poles cannot, in general, 
be used as a criterion of stability for time-varying systems (references 5 , 6 ), 
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some sufficient-existent type perturbational results exist (reference 7) which 
can provide one with some kind of feeling about the system stability behavior, 
as long as the time-varying perturbation terms are either small or slowly varying 


The unperturbed model is obtained from equation (4-2) by freezing the time 
at a particular value. If the despin torque is assumed to be applied to body 
1 and the time is frozen at t=0, then the unperturbed model is given as 


— — 


“ - 


— , — ■ 

• 

0^ 


0 0 10 



2y 




2y 

• 


0 0 0 1 



2z 




2z 

•• 





• 

0„ 


rO 

o 

o 



2y 


1 


2y 

0^ 


0 -b 0 


• 

2z 




2z 

__ _ 


— 


mm Mb 


where 

h 

^ “o 

and the time-varying perturbation comes from 




0) 

o 


which has a time constant of 0.63 t^^ for a given despin time, t^^. 


Numerical calculation shows that, for a given despin time, the time 
constant of the unperturbed system is about equal to t^^. Compared with the 
system time constant, the perturbation variation is not too slow, but is not 
too fast either. Consequently, the pole location analysis should still be 
able to reveal some information about the transient behavior of the system. 
Numerical results obtained for the unperturbed model frozen at t = 0 and t = 
tp/2 are presented in Table 4-1. The appearance of the positive poles suggests 
an unstable tendency in the time behavior. For the case when the despin 

• 

torque is applied, body 2 instability is also indicated. This is because 02^^ 
will eventually go to zero and according to equation (4-2) and the analysis 
presented in subsection 4.2.1, multiple zero poles would appear in the 0^^ or 
the 022 . process. 
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Table 4-1. OPEN-LOOP POLES/WITH DESPIN TORQUE 


tp(sec) 

POLES (t = 0) 

POLES (t = tp/2) 

20 

-0.04551 + 0.2101J 

-0.06251 + 0.1301 


+0.04551 + 0.009858j 

+0.06251 + 0. 03002 j 

60 

-0.01644 + 0.2016j 

-0.02887 + 0.1078J 


+0.01644 + 0.001 341 j 

+0.02887 + 0.007728J 

120 

-0.008304 + 0.2001J 

-0.01588 + 0.1026J 


+0.008304 + 0. 0003436 j 

+0.01588 + 0.002459J 

180 

-0.005547 + 0.2004J 

-0.01085 + 0.1013J 


+0.005547 + 0.0001 534 j 

+0.01085 + 0.001163j 


4.3 LINEAR MODEL CLOSED-LOOP BEHAVIOR 


In this subsection, the closed-loop behavior (that is, with the spring/ 
damper alignment torque) of the linear model will be investigated first with- 
out a despin torque (in subsection 4.3.1) and secondly with a despin torque 
(in subsection 4.3.2). In both cases, body 1 is assumed to have ideal 3-axis 
attitude control. Numerical calculations will be based on the nominal model 
as presented in subsection 2.3. 


4.3.1 Closed-Loop Behavior Without the Despin Torque 

For the case of the linear model with the alignment torque but without 
the despin torque, = 0 and equation (4-2) reduces to 


« 


- ■» 


— 

e 


0 0 10 


0o 

2y 




2y 

0 


0 0 0 1 


0o 

2z 

zs 



2z 

So 


-k 0 -X -b 


• 

0o 

2y 




2y 

So 


0 -k b -X 


• 

0n 

2z 




2z 

~ 


— _ 


— — 


where 


b = 


b. 

q “o 
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and 03 is a constant, 
o 

It can be shown analytically that, for any spring/damper combination 
(k > 0, X > 0) , the linearized model is always asymptotically stable, that is, 
the spring/ damper compensated system is absolutely stable which implies 
that the nonlinear model is also absolutely asymptotically stable. 


For damped harmonic motion, the critically-damped case is generally the 
most efficient for stabilization purposes. The following sample design is 
based on a simplified model for which 9^^ and 6^^ are decoupled. Choosing 


k = 0.01 = 216.64 Kg-m^/r ad-sec^ 

X = 0.2 C* = 4332.8 Kg-m^/rad-sec 

then a pole-configuration for the real system is obtained as 

= -0.1786 + 0.2275J 
P 2 = -0.17S6 - 0.2275j 
P^ = -0.0213 + 0.0271j 
P^ = -0.0213 - 0.0271j 


As a preliminary analysis, the time constant and the oscillation period for 
both the first and second conjugate pairs, (P^^, P^) and (P^j respectively, 

are given in Table 4-2 for P = - 0 + uj . 

Table 4-2. CLOSED-LOOP TIME CONSTANTS AND PERIODS/WITHOUT DESPIN TORQUE 



TIME CONSTANT 

OSCILLATION PERIOD 


T = 1/a (sec) 

T = Ztt/oi (sec) 

P-I, P 2 

== 5.6 


P 3 , P 4 

== 46.8 



From Table 4-2, it can be seen that an initial offset would damp out relatively 
fast. A more detailed analysis, along with discussion of system structure, 
zero configuration, and general design rule will be presented in subsection 4.4. 
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4.3.2 Closed-Loop Behavior with the Despin Torque 

The case of the linear model with both the alignment and the despin 
torque is represented by equation (4-2) , whereby the system coefficient matrix 
is time varying. For t > t^^, the system reduces to a model similar to that 
discussed in subsection 4.3.1, which was shown to be asymptotically stable. 
Based on a steady-state stability theorem (reference 6) then, the asymptotic 
stability of the perturbed closed-loop system with despin torque is guaranteed. 


4.4 SOME STRUCTURAL PROPERTIES, PERFORMANCE EVALUATION, AND 

GENERAL DESIGN GUIDELINE 

Owing to the existence of various combinations of two body parameters 
(mass, inertia, and joint vector), a specific set of design values based on 
a nominal model may not be very helpful. However, it is believed that 
regardless of the specific configuration of interest, the structure and some 
other properties of the compensated nominal model should be preserved. Hence, 
it is of both theoretical and practical interest, to study system structure 
from which a general design rule may be obtained. 


Basic structural properties are controllability, observability, and 

coupling/interaction. It is easy to see that the closed-loop model, with 

• • 

either ®2y’ ®2z’ ^2y ®2z ^ single output, is both (input-state) 

controllable and (state-output) observable, as long as k > 0 and X > 0 (so 

that no pole-zero cancellation can occur). In other words, each of , 0„ , 

2y 2z 

• • •' 

0- , or 0„ contain all of the four modes. 

2y 2z 


A further analysis of time-behavior requires the knowledge of zero- 
locations which in turn determines the importance of each pole. However, it 
is known that for a multi-input system, the zeros of the transfer function 
matrix are functions of both subsystem interaction and state variable feedback. 
For this underlying simple structure, it can be shown that the state output 
resolvent matrix (£(^S - A) ^ of reference 2) with 0^^ and © 2 ^ as output, has 
the form 
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®2y 


(S + A)(S^ + AS + k) + b^S 

2 A A 

S + AS + Ic 

-bk 

-bs 


^2yo 


A 

c 

A 

c 

A 

c 

A 

c 


^2yo 



A 



S“ + AS + k 


^2zo 



bk 

bs 

(S + AKs'^ + AS + k) + b S 


^2zo 


A 

A 

A 



■ 


c 

c 

c 

c 




L. 





L J 


where 

S - Laplace transform variable 

A - characteristic polynomial of the coupled model 
c 

• • 

and 0„ , » Or. , and 6„ are the initial states in the state-space. 

2 yo 2 yo 2 zo’ 2 zo ^ 


For the nominal model (see subsection 2.3), the coupling coefficient, b, 
has a value of 0.20032. Without a more sophisticated measure of system inter- 
action (reference 8 ) , it can readily be seen that the 9„ and 0 processes are 

A A z z 

weakly coupled if b is small compared with k and 1. An extreme case is given 
by b = 0 , whereby the 9^^ and 0 ^^ processes are decoupled and hence unobservable 
from each other. In addition to the disappearance of those coupling terms in 

2 A A 

the resolvent matrix, pole-zero cancellations of the roots of S + AS + k 
occur. These cancellations make each process, 0^^ or ^2z* ^°^tain only two 
modes, that is, the other two modes corresponding to say the 9 process have 
no weight in the dynamics of the 0„ process. When b ^ 0, but the subsystems 
are only weakly coupled, then the above situation should not change too much. 
This means that the roots of the characteristic polynomial of the coupled 
system are close to the roots of the characteristic polynomial of the decoupled 
system, (S'^ + AS + k) , and that only two modes are dominant in each subprocess. 


A root locus plot for the coupled system where k is held constant and A 
is varied, and where A is held fixed and k is varied is shown in Figure 4-1. 
Numerical data for this figure are given in Table 4-3. 
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Table 4-3. NUMERICAL DATA FOR FIGURE 4-1 


PLOT SYMBOL 

ic 

A 

DATA 

o 

0.01 

1.0 

-0.0096 + 0.0019J 
-0.9903 + 0.2023J 

□ 

0.01 

0.5 

-0.0174 + 0.0075J 
-0.4825 + 0.2078j 

o 

0.01 

0.2 

-0.0213 + 0.0271J 
-0.1786 + 0.2275J 

A 

0.05 

0.2 

-0.0560 + 0.1277J 
-0.1439 + 0.3281J 

0 

0.10 

0.2 

-0.0684 + 0.2176J 
-0.1315 + 0.4180J 


From the above discussion and graphic illustration, the following general 
guideline can be summarized: 

• Subprocess interaction has the effect of introducing oscillation in the 
mean time separating the two conjugate pairs. 

• Subprocess interaction is decreased as b becomes smaller, particularly 
when is small or k and X are large. 

• In general, interaction is mild because the root-locus shows a good 
command of system pole-location through spring/damper feedback. 

• The decoupled model, being a good approximation when interaction is 
weak, is suggested for use in design purposes. 

• Because of the oscillation effect introduced by coupling, slight 
overdamping is recommended. For instance, design values of k = 0.01 
and A = 0.5 may be good. However, excessive damping should be 
avoided, otherwise the slow transient conjugate pair would be dominant. 
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Section V 

A-2/A-3 CONTROL OPTION - BODY 1 
HAS X-AXIS ONLY OR NO ATTITUDE CONTROL 

The cases to be considered in this section are less restrictive than those 
discussed in Section IV. Namely, body 1 is assumed to have only x-directional 
attitude control (A-2) or no attitude control at all (A- 3) . It should be noted 
that x-axis only attitude control is implied by the condition 6^^^ = 0. The 
state variable canonical form of the linearized model to be used in this section 
is given by equation (3-9) . 

5.1 LINEAR MODEL OPEN-LOOP BEHAVIOR. DECOMPOSITION, AND INTERPRETATION 

In this subsection, the open-loop behavior (that is, without the spring/ 
damper alignment torque) of the linear model will be investigated first without 
a despin torque (in subsection 5.1.1) and secondly with a despin torque (in 
subsection 5.1.2) for both the A-2 and the A-3 control options. Numerical 
calculations will be based on the nominal model as presented in subsection 2.3. 


TR-1400 


5.1.1 Open-Loop Behavior Without the Despin Torque 

It can be readily seen from equation (3-8) that the Y and Z processes are 
unobservable from the q and ^ processes. Therefore, the state-space dimension 
can be reduced to six, and such a six-dimensional realization can be shown to 
be a minimal realization if q and ? are considered as output. Furthermore, 
if q and ? are considered as output separately, it can be shown that a five- 
dimensional realization would be sufficient and is completely observable. On 
this basis, the state-space description for the q-process for when neither the 
alignment nor the despin torque are active can be written as 



with the output 
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It should be noted that a similar form holds for the ^-process. 

For the A-2 control option, the nominal model indicates there are five 
modes in the dynamics of either n or namely 



= 0.2767j 
= -0.2767j. 


It can be shown analytically that when body 1 has only x-axis attitude control 
(6^^ = 0), there are, in general, always three zero poles, since the system 
matrix has a 3-dimensional null space. Because of the multiple zero poles and 
the observability of the system, the n or 5 process is open-loop linearly 

unstable and therefore nonlinearly unstable. It can also be shown that the n 

• • 

and ? process is not observable from the n and ? process and therefore a four- 

• • 

dimensional realization exists for the p and r, process. That is to say that 

• • 

the n and C process can be separated from the p and ? process and contains 
four modes. The imaginary poles generally imply sustained oscillation in the 

relative precession rates whereas, the multiple poles reveal, in this case, 

• • 

the nonlinear instability of the p and C process. Since the total energy is 
finite, a limit cycle type instability would be expected. 

For the A-3 control option, body 1 is assumed to be initialized with a 
perfectly matching angular velocity as compared to body 2. In so doing, the 
basic structure of the problem will remain the same as that for the A-2 control 
option. Based on the nominal model, the numerical values of the five modes 
in p or 5 become 
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P 2 = 0.2810J 
= -0.2810j 
P^ = 0. 01063 j 
P^ = -0.01063J 


Since there is only one zero pole along with two imaginary conjugate pairs, 
the linear analysis fails to predict overall stability. 


5.1.2 Open- Loop Behavior with the Oespin Torque 

The case of the linear model without the alignment torque but with a con- 
stant despin torque is discussed here relative to the A-2 and A-3 control 
modes. For x-directional attitude control of body 1 (A-2), the despin dynamics 
is similar to that discussed in subsection 4.2.2 for the A-1 control mode, 
that is, § 2 ^ should decrease linearly to zero. The same arguments as presented 
in subsection 4.2.2 are applicable here, that is, plausibility is based on 
perturbational results. Numerical results obtained for the unperturbed model 
frozen at t = 0 and tj^/2 for the case when the despin torque is applied to 
body 1 are presented in Table 5-1. Similarly, when the despin torque is 
applied to body 2, there is a complex conjugate pair of poles in the right- 
half plane for each case, which suggests an unstable tendency relative to 
misalignment . 


For the A-3 control option, body 1 is again assumed to be initialized 
with an angular velocity perfectly matching that of body 2, implying that an 
x-directional attitude control of body 1 is required to design the system. 

This reduces the problem to essentially that just discussed for the A-2 control 
mode. More generally, however, if the initial angular speed of body 1, to^, is 
different from oj^ (the value of at t = 0) , then the following despin time 
could be separated, that is. 


t^ A 


* Aw 


where A = A, A„/(A. + A_) and Aw = w - w, . 
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Table 5-1. OPEN-LOOP POLES/WITH DESPIN TORQUE 


tp(sec) 

POLES (t = 0) 

POLES (t = tp/2) 


0,0 

0,0 

30 

-0.04082 + 0.2826J 

-0.06296 + 0.1627J 


+0.04082 + 0. 005897 j 

+0.06296 + 0. 02436 j 


0,0 

0,0 

60 

-0.02104 + 0.2783j 

-0.03743 + 0.1478J 


+0.02104 + 0.001590J 

+0.03743 + 0. 009479 j 


0,0 

0,0 

90 

-0.01411 + 0.27743 

-0.02650 + 0.1432J 


+0.01411 + 0.0007181J 

+0.02650 + 0.0049013 


0,0 

0,0 

120 

-0.01061 + 0.2771J 

-0.02041 +0.14133 


+0.01061 + 0.0004062J 

+0.02041 +0.0029483 


5.2 LINEAR MODEL CLOSED-LOOP BEHAVIOR 

In this subsection, the closed-loop behavior (that is, with the spring/ 
damper alignment torque) of the linear model will be investigated first without 
a despin torque (in .subsection 5.2.1) and secondly with a despin torque (in 
subsection 5.2.2). In both cases, body 1 is assumed to have x-axis only atti- 
tude control (that is, 6^^ = 0). Numerical calculations will be based on the 
nominal model as presented in subsection 2.3. 

5.2.1 Closed-Loop Behavior Without the Despin Torque 

It can readily be shown that, without the despin torque and using the 
alignment torque as defined by equation (3-4) , the state variable canonical 
form of the linearized model can be expressed as 


• 

n 


0 

0 

0 

0 

1 

0 


n 

? 


0 

0 

0 

0 

0 

1 



Y 

= 

0 

0 

0 

-Pi 

0 

P 2 


Y 

Z 


0 

0 

Pi 

0 

-P 2 

0 


Z 

• 

n 


-k 

0 

0 

P 3 

-\ 

-P 4 


n 

c 


0 

-ic 

-P 3 

0 

P 4 

-X 


f — 
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where k = hk and X = hA . This system is composed of two coupled subsystems, 
namely the ri and 5 process and the Y and Z process. Only the n and C process 
can be directly controlled by the spring/damper alignment torque. A further 
study of the system reveals that the Y and Z subsystem has associated with it 
two pure imaginary poles which will introduce some oscillation throughout the 
system. Therefore, a pair of conjugate poles are expected to exist which can- 
not be moved very far from the imaginary axis. However, it is not unreasonable 
to assume that these two slow-transient poles will have little effect on n or 5 
if the coupling is weak. An analytic proof of the absolute stability of this 
linear model would probably be difficult to obtain because of the system order. 
In Section VI, a Liapunov-type approach involving the original nonlinear system 
will be presented. 

In the following sample design, the spring 
to make the isolated n and ^ process critically 

k = k/h = 0.01/h = 85.601 kg-m"^ /rad-sec^ 

X = x/h = .2/h = 1712.0 kg-m^/rad-sec 

results in the following pole-configuration: 

= -0.1718 + 0.29l7j 
= -0.1718 - 0.2917j 
= -0.02604 + 0.03537j 
P^ = -0.02604 - 0.03537j 
= -0.002088 + 0.02032j 
Pg = -0.002088 - 0.02032j 

It is readily recognized that the first two pairs are attributed primarily to 

• • 

the ri and ^ subsystem and the third pair is due primarily to the Y and Z sub- 
system. The first pair represents fast- transient modes, whereas the second 
pair is slower, having a time constant of 38.4 seconds and a modal oscillation 
period of 180 seconds, which implies a pretty good damping effect. The third 
pair corresponds to very slow transient modes which have a time constant of 
478.7 seconds and an oscillation period of 390.3 seconds. 


and damper constants are chosen 
damped. Choosing k and X as 
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Although it is reasonable to expect the third pair to have little impor- 
tance in n/n > Ti/C , Ti/n , n/? , etc. , it may have some persistent effect in 

• • 

the impulse reponse of etc. A detailed analysis would require a 

pole-zero map and the understanding of the system structure, both of which will 
be presented in subsection 5.3. 

It might also be expected that the effect of the slow transient pair on 
etc. could be reduced by an increase in the energy storage and dissipa- 
tion capability, that is, by increasing k and X. Choosing k = 1 and X = 2 
gives k = 8360.1 and X = 17120, resulting in the following pole-configuration; 

= -1.33505 + 0.49056j 

= -1.33505 - 0.49056j 

P^ = -0.664929 + 0.243993j 

P^ = -0.664929 - 0.243993j 

P^ = -0.0000133094 + 0.030170j 

P, = -0.0000133094 - Q.030170j 
0 

It is Interesting to note that the third pair moves even closer to the imaginary 
axis. This phenomenon will be self-evident after the discussion of interaction 
in subsection 5.3. From a pole- con figuration standpoint, there does not appear 
to be any improvement relative to an increase in k and X. However, a more 
detailed pole-zero and time response study to be presented in subsection 5.3 
wil3. show that the latter of the above two designs is indeed much better than 
the first. 

5.2.2 Closed-Loop Behavior With the Despin Torque 

The case of the linear model with both the alignment and the despin torque 
is represented by equation (3-9) whereby the system coefficient matrix is time 
varying. For t > t^, the system reduces to a model similar to that discussed 
in subsection 5.2.1, whereby asymptotic stability of the perturbed closed-loop 
system with despin would be preserved (based on a steady-state stability theorem 
(reference 6)). 
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5.3 SOME STRUCTURAL PROPERTIES, PERFORMANCE EVALUATION, AND GENERAL 
DESIGN GUIDELINE 

It can be readily shovra that the spring/damper compensated system for any 

k > 0 and X > 0 is completely observable with the six-dimensional state-space 

• • 

representation and the scalar output of n> n or 5 . Therefore, all of the 
six modes will appear in ri or C- It was shown in subsection 5.2 that there is 
one conjugate pair corresponding to slow-transient modes and for which a root- 
locus analysis showed that they cannot be moved far away from the imaginary 

axis. Since we are interested in the dynamics of n and ^ as opposed to that 

• • 

of Y or Z, the zero locations and hence the residues of those modes in the 
resolvent matrix of n and 5 are essential in determining the importance of 
those slow-transient modes in the dynamics of n and 5 . Prior to giving any 
numerical results, a general discussion of the system interaction will be 
presented, which should be useful when considering the variant models. 

It was mentioned in subsection 5.2 that the overall system can be viewed 
as consisting of two subsystems coupled together where only one of these sub- 
systems could have external inputs (see Figure 5-1) . Mesarovic (reference 8) 
labels this form as a V-canonical structure. 



Figure 5-1. V-CANONICAL STRUCTURE MODEL 
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These two subsystems can be denoted by 


and 


[H^(S)] 




"1 


where 


[H^(S)] 


• • 

resolvent matrix from the n , ? > n and C sub space to n and C* 

o o o o 

« • 

resolvent matrix from the Y and Z subspace to p and 

o o 


It is interesting to note that the subsystems are coupled through the parameters 

• • 

p- and p„. If p- is zero, then the Y and Z processes would be unobservable from 
^ ^ • • 
the n and ^ processes and hence the two open-loop poles of Y and Z would be 

cancelled by two zeros in [H^(S)]. In general, there are two zeros in [H^(S) ] 

which are zeros having the generic form of 

(S^ + p^) -h f(p3,S) 

where f(p 3 ,S) is a rational function In the Laplace Transform variable S which 

is continuous in p^ and equals zero if p^ = 0. Also, the two open-loop poles 

of Y and Z are given by +p^^j and -p^^j . If the interaction is weak or p^ is 

small, there will be two zeros in [H (S)] which are close to these two slow- 

transient poles. Consequently, they will have little influence on [H^(S)] and 

the behavior of [H (S) ] should be similar to that discussed in subsection 4.3.1 

o' 

for the closed-loop system without the despin torque and body 1 with 3-axis 
ideal attitude control. 


as 


On the other hand, if P 2 0» the resolvent matrix [H^(S)] can be x<7ritten 


IH^(S)] = 


S" + AS + k "^4^ 


— “ 
^1 ^3 ^3^ 

A A 


o2 . 2 r.2 . 2 

0 0 


S + p^ S + p^ 

^4^ + AS + k 


P3S -P3P3 

A ^ 


^2 . 2 n2 . 2 

0 0 


S + S + P2^ 
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P3[pi(S^ + AS + k) - p^S^] 
A^(S^ + pp 

p^S^^S^ + AS + k + Pj^ p^j 
A (S^ + p^) 

O 1 


PgS[s^ + AS + k + pj 

~P3[pj^(sV AS + k) - p^S^] 
_ _ 

AqCs^ + pp 


where A is the open-loop characteristic polynomial of the n and z processes. 
The structure of this configuration is expressed as shown in Figure 5-2. 



Figure 5-2. SIMPLIFIED STRUCTURE I 


Considering the six poles according to their order as discussed in sub- 
section 5.2.1 (that is, the fast-transient pair as the first pair, the 
intermediate-transient pair as the second pair, and the slow- transient pair as 
the third pair) , then the transformation for p^ = 0 indicates that both the 
second and third pairs will be important. Even when P 2 ^ Q, the situation 
should not change too much if the interaction is weak or intermediate. Also, 
as P 2 gets smaller, the slow- transient pair will move closer to the imaginary 
axis and as p^ gets smaller, the slow- transient pair will be less involved in 
the n and Z dynamics. These observations are summarized in Table 5-2 where 
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the important pole pairs were determined by the residue of the respective 
transfer function. Analogous conclusions can also be made for the dynamics 
of 


Table 5 - 2 . IMPORTANT POLE-PAIRS 


TRANSFER FUNCTION 

IMPORTANT POLE-PAIRS 

n/n. 


2nd 

pai 

r 

0 

n/4o 


2nd 

pa1 

r 

n/n 


and 

2nd 

pai r 

U 

h/Cq 

^st 

and 

2nd 

pair 

n/Y 

2nd 

and 


pair 

0 

2nd 

and 

3rd 

pair 


Several points relative to the above discussion should be mentioned at 

this time. First, the above conclusions are more plausible when interaction 

becomes weaker (when p^ becomes smaller or k and X become larger as compared 

with p_) . Secondly, the slow-transient modes are Important only in the n/Y 
. . . . . . ® 

and n/Z transfers where Y and Z are given relative to 0. , 0_ , 0. 

o o o ° lyo 2yo Izo 

and ^2zo assumed fairly small. Thirdly, the important pole-pairs 

listed for a particular transfer function were determined through a comparison 
with that transfer function only and the Importance of each pole-pair relative 
to all the transfer functions listed in Table 5-1 was not considered. Finally, 
the amplitudes for each pole-pair in every impulse response were not determined. 

a 

However, it is expected that the amplitudes of the impulse responses n/Y and 

. o 

n/Z would be reduced when interaction becomes smaller, 
o 

The coupling factors Pj^, P 2 * ^4 nominal model when body 1 

has ideal x-dlrectlonal attitude control (the A-2 control mode) are numerically 

• • 

given in Table 5-3. It is interesting to note that the Z, Y to n, 4 coupling 

factor Pj la pretty large if k and X are chosen as ic - 0.1 and X - 0.2, 

• • 

so that a larger influence from Z^, Y^ to n, 4 is expected. In Table 5-4, 
numerical data is given for the poles, zeros, and gain associated with the 
appropriate element of the resolvent matrix for each of the Indicated transfer 

functions. In each case, the spring and damping were chosen as k ■ 0.01 and 

- -4 

X * 0.2. In addition, any number having an absolute value less than 10 was 
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truncated. It is of interest to note that these numerical results agree quite 
well with the summary of general characteristics. 


Table 5-3. COUPLING FACTORS 



COUPLING 

FACTORS 



Y to Z 

= 0.03039 



C,n to 'J'.Z 

P 2 = -0.02118 



Y,Z to s,n 

p^ = -0.03534 



n to c 

P 4 = 0.02463 



Table 5-4. POLES, ZEROS, 

AND GAIN OF THE RESOLVENT MATRIX 

TRANSFER 

POLES 

ZEROS 

GAIN 

n 

% 


-0.1718 + 0.2917J 

-0.1772 + 0.2792J 



-0.02604 + 0.03537j 

-0.02924 

1 


-0.002088 + 0.02032J 

-0.008110 + 0. 02261 j 


n 

-0.1718 + 0.2917J* 



^0 

-0.02604 + 0.03537j 


0.002463 


-0.002088 + 0.02032J 

0,0 • 


n 

-0.1718 + 0.2917J 



■^0 


-0.02604 t 0.03537J 

-0,09769 + 0.08614J 

1 


-0.002088 + 0. 02032 j 

-0.002309 + 0.02322j 


n 

-0.1718 + 0.2917j 



• ' 

r 

^0 

-0.02604 + 0. 03537 j 

0 

-0.2463 


-0.002088 + 0.02032j 

0,0 


n 

-0.1718 + 0.2917J 



Yo 

-0.02604 + 0.03537j 

-0.01098 + 0.031 26 j 

-0.09780 


-0.002088 + 0. 02032 j 



n 

-0.1718 + 0.2917J 

- 0.1 

' 

Zo 

-0.02604 + 0.03537j 

- 0.1 

-0.3534 


-0.002088 + 0. 02032 j 

0 
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In general, the time response of a given transfer function can be expressed 


as 


r(t) = G 


y R. e x+2 y C. ex cos(w.t + 6,) 

X X X X 

x=l x=l 


m 


-X.t 


where 

- the residue corresponding to the real mode, 

2C^ - the residue corresponding to the complex modes, 

6^ - the phase shift 

G - the gain 

In Table 5-5 numerical data are given for k = 0.01 and X = 0.2 summarizing the 
residues and phase shifts for each of the modes considered. 


Table 5-5. RESIDUES AND PHASE SHIFTS 


TRANSFER 

POLES (PAIRS) 

2GC^- (WEIGHT) 

•PHASE SHIFT 

n/n^ 

^st/jnd/jrd 

0.08055/0.74146/0.25925 

-0.9288/-0.2455/-0.4596 



0.08043/0.74041/0.25898 

-5.6412/-1 .8T53/-1 .1111 

0 

n/n„ 

l=V2"“/3'"‘' 

2.7276/3.257/0.5297 

-1.9673/-1.1806/-1.9279 

0 

0/^0 

^st^jnd/srd 

-2.7276/-3.257/-0.52968 

-3.5331/0.39014/2.7843 



-3.0814/-16.2723/-18.2039 

-5.6905/-1.5577/-1.7234 

vz: 


-3.081 4/-16.2719/-18.2036 

-4.1197/-3.1285/-0.15268 


Numerical calculation shows that the time impulse response of jp/Z^I has 

maximum peak values of 12.4, 13.23, 9.44, and 6.83 radians at times of 40, 155, 

310, and 465 seconds respectively. Similar peak values also appear in the 

|n/Y I time response. Consequently, initial wobble-may have a large effect on 

the relative misalignment. The nominal value of 0.1 rad/sec for 6„ is 

2x0 

almost the largest possible spin-rate. For a mild wobble, a value of 0.01 
rad/sec for 6^ ^ ®2zo reasonable, resulting in a damped sinusoidal 

oscillation In misalignment with an amplitude of about 3.159 degrees during 
the first 200 seconds. This oscillation, having an initial period of about 
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250 seconds which would increase with time, has its amplitude decreasing 
exponentially with a time constant of about 500 seconds. Although the 

situation for this kind of initial wobble is not too serious, an undesired 

♦ • 

situation may occur if either ^2yo ^2zo ^ value of 0.05 rad/sec since 
this would result in an initial modal oscillation amplitude of 15.79 degrees. 

Some attention then should be paid to the behavior of the r)/Y and n/Z 

o o 

transfers. It has been observed that the two slow- transient modes are unavoid- 
able regardless of the values chosen for the spring and damper coefficients. 

So at best, one can only try to reduce the amplitudes of n/Y^ and ^/Z^. At 
this preliminary design stage, it is desirable to develop a formula which can 

be used to estimate the oscillation amplitude of n/Y or n/Z for a given 
_ o o 

value of k and X. Based on the simplified model presented below, such a for- 
mula is developed, the accuracy of which improves as the subsystem interaction 
becomes weaker. 

Suppose that the n and 5 processes were decoupled and would 
• • • • 

the Y and Z processes, but the Y and Z processes could influence 

processes. Graphically, this situation is illustrated in Figure 

• • • • 
to the structure of the Y and Z subsystem, the resulting Y and Z 

are sinusoidal functions with a frequency of p^/2ir, that is, 

Y(t) = Y cos p- 1 + Z sin p,t 
o *^1 o *^1 

Z(t) = -Y sin p, t + Z cos p.t. 

o *^1 o *^1 

Assuming that the spring/damper feedback loops are closed for the n and C 
processes, the above problem is readily recognized as a problem of forced 
oscillation with damping, which can be described by 

X + Xx + kx = a cos cot 

and the particular (steady state) solution to this equation can be written as 
X 2 (t) = A cos (cot - (()) 

where the amplitude A and phase <j) are given by 


not influence 
the n and z, 
5-3. Due 
pro cesses 
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A = 


^ 2.2 , .2 .2 


0) ) + X ta) 


and 



Consequently, it is readily seen that for Y = 0, an estimate of the amplitude 

o • 

of the slow-transient modes in the impulse response n/Z^ can be given by 

IPqI 

A = ; i ■_ (5-1) 

2,2 ^ 2 

V(k - p^) + A p^ 

From this formula, it is easy to see how the amplitude will be reduced as k 
and A increase, a phenomenon which has been predicted from a subsystem 
interaction point of view. 


■-Ay 
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Figure 5-3. SIMPLIFIED STRUCTURE II 
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A comparison of the estimated amplitude and the calculated "true" amplitude 

is presented in Table 5-6. From this table it is immediately apparent that as 

k and A increase, the accuracy of the estimate provided by the formula increases. 

In addition, the data indicates that by using more control effort (k and \) , the 

overall system performance is maintained while reducing the amplitude of the 

p/Z and n/Y transfers. 

'' o o 


Table 5-6. A COMPARISON OF ESTIMATED AND ACTUAL AMPLITUDES 


CONTROL 

EFFORT 

ESTIMATED 
AMPLITUDE. A 

TRUE AMPLITUDE FOR THE 
3 rd PAIR IN n/i^ 

TIME CONSTANT FOR 
THE 3>"a PAIR (SEC) 

OVERALL AMPLITUDE FOR 
IMPULSE RESPONSE n/Z^ 

R = 0.01 
1 = 0.2 

32.35 

18.20 

478.76 

13.23 

R = 0.36 
1 = 1.2 

0.9791 

0.9414 

17361 

0.9348 

R = 0.64 
X = 1 .6 

0.5514 

0.5391 

39308 

0.5367 

R = 1 
1=2 

0.3530 

0.3463 

75187 

0.3450 


The root locus diagram shown in Figure 5-4 shows the general trend of the 
system-pole movement. It should be noted that the third pair of poles cannot 
be moved further to the left from the imaginary axis by increasing the damping 
ratio. Although a slight movement to the left was observed for smaller values 

of k and A, decreasing the values of k and A would also result in undesirably 

• • 

larger amplitudes of the slow-transient modes in the p/Z or rj/Y transfer. In 

o o 

general then, one should not pay much attention to the location of these two 
poles. Instead, the formula should be used to choose suitable magnitudes of 
Ic and A which will give satisfactory performance for the first and the second 
pairs and yet make the amplitude of the third pair in the h/Z^ and h/Y^ transfers 
smaller than a maximum tolerable value. 
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Figure 5-4. ROOT LOCUS DIAGRAM 


At this point, one can state some generalizations and a design guideline 

for the spring/damper compensated coupled spinning system: 

• • 

1. The influence of the Y, Z subsystem on the p and C processes cannot 
be efficiently controlled by spring/damper feedback. In particular, 
two slow-transient modes which are primarily attributed to the Y, Z 
subsystem will have considerable importance in the h/Zq, p/Zqs 

and c/Zo transfers. The importance of these slow-transient modes 
will be reduced as the Y, Z to p, C interaction becomes weaker. There 
are two ways to reduce the interaction. One is to reduce P 3 bilinearly 
by adjusting the body 1 x-directional spin-rate (to be discussed in 
subsection 5.4). Another one is to increase the magnitudes of k and A. 

A formula which gives an approximate relation between k, A, and the 
amplitude of the slow-transient modes was given by equation (5-1) . 

• • 

2. The coupling of the Y, Z and p, C processes has the effect of introducing 
some oscillation in the mean time and tends to separate the first pair 
from the second. 

3. Based on items 1 and 2, a stabilization design procedure can be out- 
lined in the following three steps: 

a) Use of equation (5-1) to choose the magnitudes of k and A which ^ 
would meet a maximum tolerable limit on the amplitude of the p/Z 
transfer. 
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b) Use a decoupled/non-interacting model to choose^ the damping ratio, 
etc., to control the behavior of n/ho, h/Co> h/no» u/So» C/ho> 

and C/Co* Critical or slightly over damped is recommended, 

c) Run a root-locus and time-response analysis, etc., to confirm that 
this preliminary design is satisfactory. 


5.4 EFFECTS OF INITIAL ANGULAR SPEED MATCHING 


In a number of instances in the foregoing discussion, ideal x-directional 
attitude control on body 1 was assumed. If in addition, pre-spin docking is 
assumed, then some effects on the misalignment dynamics might be expected, 
because the strength of system interaction could be changed. Imposing the 


initial condition of perfect angular speed matching, that is, = ®2xo 

on the definitions of the coupling factors p^ and p^ given in subsection 3.3 




as 


P 

P 


2 

3 




(a2 - 3 ^ Aj 


results in 

P2 “ “o‘®2 \ ■ ““l “2^^ 

p^ = w^(a2 - a^ 


If p^ and p^ are denoted by p^ and p^ for the case when body 1 has x-directional 

A 

attitude control and and p^ for when body 1 has no attitude control, then 
for the initial condition of perfect angular speed matching, the conditions 




hold as long as a 2 A^^ < 2a^A2. The implication here is that the n, ? and 

• • 

Y, Z interaction would be weakened, thereby providing better overall system 
performance. Based on the nominal model, a 2 A^^ and a^A 2 are equal to 364,139,780.7 
and 4,319,741,524 respectively and the above inequality is seen to be satisfied. 

The above discussion suggests that for any given body 2 (satellite) spin 
rate (w^) , there always exists a body 1 (tug) matching speed (w*) given by 
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(l) 

o 


such that the n» ? and Y, Z processes become noninteracting. Therefore, the 
two slow-transient modes will not affect the dynamics of n and c;. It is 
interesting to note that with this ideal matching speed, the model reduces to 
a form similar to that found for the case where body 1 has ideal 3-axis attitude 
control. Therefore, the same performance of relative misalignment as for 
the case when body 1 has full attitude control can be achieved by prespinning 
body 1 to the ideal matching speed. However, this kind of bilinear control is 
not of practical interest, at least for this nominal model, because of the 
large value of the ratio, a.^ky 
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Section VI 

ABSOLUTE STABILITY AND THE STABILITY REGION OF THE 
SPRING/DAMPER COMPENSATED COUPLED TWO-BODY SYSTEM 

As was mentioned in Section II, the nature of the coupled two-body problem 
is nonlinear even though it was approached by formulating and investigating a 
linear model. This linearization approach enables one to take full advantage 
of the simplicity of linear structure for system design and was justified on 
the basis of Liapunov's first method, that is, linear asymptotic stability im- 
plies nonlinear asymptotic stability. However, it should be emphasized that 
some properties, such as the global property of linear stability generally may 
not be carried over to the nonlinear model. In fact, nonlinear asymptotic 
stability is, in general, a local property, that is, not all initial states will 
result in an asymptotic stable solution. The maximum region in the state space 
for which every initial state would end up with an asymptotic stable solution 
is called the (asymptotic) "stability region". Estimation of the stability 
region is a very important and practical engineering problem even though there 
are no classical stability theorems which can solve this problem. Liapunov's 
second method, however, can provide some information, but the "estimation" ob- 
tained is strongly dependent on how good a choice is made for the Liapunov 
function. Generally speaking, the estimated stability region is conservative. 
One systematic approach (reference 2) is to utilize the quadratic Liapunov 
function (V) based on the linear model. It can be obtained by solving a 
Liapunov equation and then finding its time derivative (V^^) with respect to the 
nonlinear model. An estimation of the stability region is then obtained by 
finding the largest region in the state space for which V is positive and is 
negative (see Figure 6-1) . This search-type scheme can be formulated as a 
mathematical programming problem. The estimate of stability region based on V 
is given as 

P. = {x I V(x) < K*} 

where 

and 

» = (x) = 0} . 
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As the nonlinear terms become more complicated such as with the problem herein, 
the above procedure becomes increasingly more difficult. Sometimes, however, 
the nonlinear model could be globally asymptotic-stable, such as the one- 
dimensional Lifenard-type damped harmonic motion with a nonlinear spring 
(reference 6) which can usually be revealed by physical reasoning. In the 
following, a heuristic proof on the absolute stability and global stability 
will be given, which uses physical and system-theoretic reasoning. 





Figure 6-1. STABILITY REGION ESTIMATION BY QUADRATIC LIAPUNOV FUNCTION 

First, it is easy to see that all the solutions are bounded. This is 
because the total energy is finite. Secondly, it should be noted that it is 
impossible for ^n equilibrium point with n 0 and 5 0 to exist. This will 

become apparent from the following discussion. First, consider the system 
having the generic form as indicated in Figure 6-2 where the interacting 
channels may contain some kind of nonlinearity. In mathematical form, one can 
write 

n = -kn - Xn + f^(n, n, 4 , y, z) 

? = -kc - + f^(n, n, c, c, T, z) 
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Y = g]^(n, n, 4, i, z) 

Z ■" 82^^* Cj » z) • 


Now, when n = n and ? = 4 where both n and 5 are nonzero constants, then 
- 00 00 

n = 0 and r, = 0. This implies that 6„ = - 9 and 0. =9, =6 and 

^ 2y ly y 2z Iz z 

consequently Y ~ 6 and Z = 0 . There are two possibilities for Y and Z: 

one is that Y and Z are constants and the other is that Y and Z are time 

varying. An extreme case of the first category will be examined first, 

• • • • • • 

namely Y = 0 and Z = 0 which implies that 0„ =9, =0 and 0^ =9-, =0. 

2y ly 2z Iz 

Physically, this is a situation for which the relative misalignment between 
the two bodies is not zero and in addition each body is not rotating about the 
body y and z axes. Due to the existence of the spring/damper compensation, 
it is very plausible that this cannot be an equilibrium situation in the 
p-? plane. 



Figure 6-2. GENERAL NONLINEAR STRUCTURE MODEL 
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• • 

Now suppose that Y and Z are equal to nonzero constants such that 


”1 


~ " 

n 


n 

_ 

o 



4 



r 

o 


• • 

along with Y and Z form an equilibrium state. Accordingly, a contradiction 
o o 

to the case with Y = 0 and Z = 0 is followed by Galilean invariance (reference 

9). Next, consider the situation where Y(t) = Y*(t) and Z(t) = Z*(t) , that is, 
• • • • 

Y and. Z are not constants. In order to satisfy n = 0 and ? = 0, then one 

must have both 

fiCnQ. 0, 0, Y*, z*) = k 

f (n^, 0, t , 0, Y*, Z*) = k s 
Zoo o 


and 


Y* = 0, 0, Y*^ 

z* = §2^%’ ^o’ 


Note that under this assumption, the solution to the first two of these algebraic 
equations should not be unique. In particular, some constant solutions for 
Y and Z do exist. For example, if Y = Y and Z = Z are taken as constant 
solutions, then 


Hi 

!-■ 

o 

0, 


• + 

y', 

• + 

z^) 


o 

o 

and 



•t 

Y^, 

. + 

Z^) 




0, 

c„, 0, 


0, 


0, Y*, Z*) = kq 

o 

0, Y*, Z*) = k? . 

o 


In other words, one can imagine that the n, K subsystem has constant inputs 
as indicated in Figure 6-3. However, this turns out to be the same situation 
as the first category which has already been shown to be impossible. So, the 
nonexistence of equilibrium points with nonzero misalignments has thus been 
shown . 
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Figure 6-3. A DECOMPOSED VIEW 


Finally, the possibility of a limit cycle type instability is considered. 
Suppose that the stability is not global and that an initial state outside the 
stability region might result in a trajectory oscillating along a limit cycle. 
Owing to the energy dissipation through damping, the oscillation period must 
become longer and longer with both |ri| and |cj going to zero as time goes 
to infinity. This situation then approaches the situation discussed before, 
that is, the nonzero misalignment equilibrium point which has been shown to be 
impossible. 

In the above discussion, it has been shown that any initial state would 
result in zero misalignment in the q , ^ space, that is, the stability is global. 
Meanwhile, no particular values of k and X were assumed so absolute stability 
is also proven as a general result. To look at absolute stability for the case 
where body 1 has x-directional attitude control, very small spring/damper con- 
stants were selected. These numerical results are presented in Table 6-1 from 
which it can be seen that, as long as k > 0 and X > 0, then the system is (at 
least in principle) asymptotically stable. 
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Table 6-1 . POLES FOR SMALL SPRING/DAMPER CONSTANTS 


SPRING/ DAMPER 

POLES 

k = 0 

A - 2. X 10"^ 

0,0 

-0.17803 X 10"^ + 0.27674j 
-0.21965 X 10'^ + 0.14030 x 10"^j 

k = 1. X 10’^^ 
A = 2. X 10"^ 

-0.46185, - 0.56772 
-0.17803 X 10"^ + 0.27674J 
-0.21108 X 10“^ - 0.21993J 

l< = 1. X 10“^° 
A = 2. X 10“^ 

-0.10981 X 10"^ + 0.31266 x 10"^j 
-0.17803 X 10“^ + 0.27674J 
-0.10984 X 10“^ + 0.31268 x 10"^j 

i< = 1 . X 10“^° 
A = 2. X 10’^ 

-0.10976 X 10"^ + 0.33138 x 10"^j 
-0.17803 X 10"^ + 0.27674J 
-0.10984 X 10"^ + 0.33141 x lO'^j 

o 

I — 

o 

X 

r— O 
It II 

1^ l'< 

0 + 0.33138 X 10"^j 
0 + 0.27674J 
0 + 0.33141 X lO'^j 
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Section Vi I 

SIMULATION RESPONSE STUDY 

7.1 GENERAL 

A digital computer simulation called 2B0DY (reference 10) was developed 
for studying the post-docking behavior of the Space Tug during recovery of a 
spinning satellite. The simulation accounts for the translational as well as 
the rotational motion of both the Tug and a target vehicle. The two vehicles 
are modeled as rigid bodies physically attached at a single point about which 
relative rotation between these two bodies is permitted. In addition, the 
target vehicle is allowed the freedom to spin about an axis in the proximity 
of the docking port axis. Such a configuration is analogous to the soft-dock 
regime wherein capture or initial latch has been achieved but final latch or 
hard-dock has not. Options relative to vehicle dynamics include ideal 3-axis 
attitude control of the chase vehicle, a spring/damper type alignment torque 
applied at the joint between the vehicles, and a despin torquer for despinning 
a spinning payload. These nonlinear models are described in more detail in 
Appendix D. This simulation, providing both digital and graphical output was 
used to study the dynamical behavior of various post-docked configurations as 
well as to verify the results of the analytical stability presented previously 
in Part I, Sections II through VI. 

7.2 SIMULATION RESULTS 

Simulation studies were performed to answer several questions concerning 
post-dock behavior during the recovery of a spinning satellite. The primary 
concern was to determine whether the post-docked configuration is stable and 
well-behaved for various: 

• target inertia characteristics 

• alignment spring/damper characteristics 

• despin torquer limits. 

Figures 7-1 and 7-2 show the overall results of these studies for an 
uncontrolled Tug and for one with perfect attitude control respectively. 

Three target inertias are considered: a disk* an end-dock cylinder and a 
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Figure 7-2. TWO-BODY BEHAVIOR FOR TUG 3-AXIS ATTITUDE CONTROL 
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side-dock cylinder. In each case the target is assumed to be initially at a 
constant angular rate about the docking port centerline with some misalignment 
and wobble. The response curves show qualitatively the parametric angular 
motion of the target centerline in inertial space and the time response of the 
angle 3 between the Tug and target vehicle centerlines, with and without 
alignment springs and dampers . 

The results in Figure 7-1 show that the uncontrolled post-dock configu- 
ration is stable for all three types of targets, provided alignment springs 
and dampers are present. The motion of both centerlines is well behaved and 
the alignment angle tends toward zero in all cases. The righthand column in 
Figure 7-1 illustrates the fact that the post-dock configuration tends to be 
unstable without alignment springs and dampers. The angle between the two 
centerlines diverges indicating the tendency of the two vehicles to pitch and 
yaw into each other. Figure 7-2 shows that alignment springs and dampers are 
still necessary even when the Tug has ideal attitude control. 

Alignment damping and natural frequency were parameterized as was despin 
torque to determine appropriate settings. Figure 7-3 shows the behavior of 
the misalignment angle 3 for several damping ratios with all other conditions 
held constant. The conclusion is that critical damping (? = 1.0) is best 
in driving misalignment quickly to zero with a minimum of oscillation. Natural 
frequency f^ should be as high as possible for fast settling time, as shown in 
Figure 7-4, but will be limited by the weight of the springs in relation to 
vehicle weights . 

The rate of despin is linearly proportional to despin torque. A despin 
torque of 100 ft-lb can despin the target from 0.1 rad/sec to zero in approxi- 
mately 43 seconds and one third of a revolution. Lower torques require longer 
times to despin as shown in Figure 7-5, but in no case is there a stability 
problem, provided alignment springs and dampers are used. Despin torque was 
applied about the centerline of the Tug as well as the target vehicle with 
little difference noted in the response. The conclusion is that the exact 
direction and magnitude of despin torque are not critical. 
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BODY 2 X-AXIS COMPONENT OF ANGULAR VELOCITY (RAD/SEC) 



Figure 7-5. DFSPIN TORQUE EFFECT 
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» An alignment torque damping factor in the proximity of critical 
damping Is best for reducing misalignment between the docking port 
asses to Kero with a minimum of oscillation. 

i An alignment torque natural frequency should be as high as possible 
for fast settling time, but will obviously be limited by structural 
weight eonwlderatious. 

• The despill time is linearly proportional to the despin torque (for 
a reasonable range of values) with a smaller torque requiring a 
longer desplii time, 
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Table 8-1. CONFIGURATIONS STUDIED 


TUG/PAYLOAD 

CONFIGURATION 

ANALYTIC 

STUDY 

NUMERIC 

STUDY 


B 

DISK 

PAYLOAD 


/ 




CcrT<^ CTD 

END-DOCK 

PAYLOAD 


6 

X 

SIC 

p/ 

)E-DOCK 

lYLOAD 


*/ 

- 
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Table 8-2. ANALYTICAL STABILITY STUDY RESULTS 


ALIGNMENT 

TORQUE 

DESPIN 

TORQUE 

TUG ATTITUDE CONTROL 

NONE 

X-AXIS ONLY 

3-AXIS 

NO 

NO 

NO CONCLUSION (FOR 
NONLINEAR MODEL) 

UNSTABLE 

NO CONCLUSION (FOR 
NONLINEAR MODEL) 

YES 

NO 

ABSOLUTELY STABLE 

ABSOLUTELY 

STABLE 

ABSOLUTELY STABLE 

NO 

YES 

UNSTABLE 

UNSTABLE 

UNSTABLE 

YES 

YES 

ASYMPTOTICALLY 

STABLE 

ASYMPTOTICALLY 

STABLE 

ASYMPTOTICALLY 

STABLE 


.S' 

o 

o 
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It may be concluded that, for the cases studied, the two-body coupled 
system is stable as long as some dissipative/energy storage (spring/damper) 
alignment torque exist. In addition, no system instability is encountered 
for reasonable values of despin torque. 

In the foregoing studies, it was found that the results derived from the 
nonlinear simulation agreed well with those predicted by the linear analytical 
stability analysis. However, it should be understood that the linear analytical 
model and the nonlinear simulation model have different functions. And because 
of their inherent advantages and disadvantages, a judiciously combined utiliza- 
tion of the two provides the greatest amount of informat ’’on and is therefore 
essential to an efficient and successful analysis and design. 

In general, the analytic model has the advantages of providing not only 
stability Information (except for a few open-loop cases) , but also a design 
guideline, an understanding of system structure, and therefore possibly 
parametric effects. The disadvantages include a loss of generality (for 
instance, there is no simple exact linearized model for the side-dock payload 
configuration depicted in Table 8-1) and a loss of time fidelity (nonlinearities 
are absent). On the other hand, the simulation model has the advantages of 
providing both more genuine results (nonlinear effects are taken into account) 
and information for more complicated cases (like the above mentioned case 
involving a side-dock target). Also, it is suitable for parametric studies. 
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Appendix A 

GENERAL EQUATIONS OF MOTION FOR TWO COUPLED RIGID BODIES 


The general rotational equations of motion for two rigid bodies connected 
at a single point as shown in Figure D-1 is given by a (6 x 1) matrix equation 
as 

J w = A 


where 


J 


and 


1 





IJ , K 

^ 

II 

• 3 


. = 


' I* 
1 2 


(I) 

2 


*2 

mt 1 — < 


- - 


— 


I* = I. - 

1 X 


y d . d . , i=l , 2 

XX 


K = y dj^ T' d^ 


^1 = ^1 '^l ^1 " “l ^1 ‘"l 


^2 ^2 ^2 “ ^2 ^2 ^^2 ^2 “2 


F* 

1 


m 

y 


ra. 


m„ 

_J. 

T' f_ 

2 

- — f 

m 

2 

m 


+ m_ 


1 

2 


™1 

m^ / 



2 "2 2 


y 0,^ d^ 


- (3x3) inertia matrix for body i, i = 1, 2 

d^ - (3x1) column matrix of the joint vector components for body i, 

i ~ 1> 2 (the joint vector is from the center of mass of body i to 
the hinge point between the two bodies) 


T* - (3x3) coordinate transformation matrix resolving a vector defined in 



“ (3x1) column matrix of the angular velocity components of body i, 
i = 1, 2 

“ (3x1) column matrix of the total external torque components acting 
on body i, i = 1, 2 
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f. - (3x1) column matrix of the total external force components acting on 
^ body i, i = 1, 2 

L. - (3x1) column matrix of the joint torque components acting on body i, 
^ i = 1, 2 

It should be noted that the tilde (~) placed over the symbol for a column 
vector indicates the 3-dimensional vector cross product operator, which for 
any vector 



is defined as 



A- 2 
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Appendix B 

SIMPLIFIED LINEARIZED EQUATIONS OF MOTION 
FOR TWO COUPLED RIGID BODIES 


In Appendix A, a summary of the completely general nonlinear equations 
of motion for two coupled rigid bodies was given. In the appendix herein, 
these equations of motion will be simplified according to the following 
assumptions : 

• No external forces or torques are assumed to act on either of the 
two bodies. Therefore, 


f s= f = 0 * 1= 0 * = 0 

ff ^1 ^2 


and the expressions for and F* in Appendix A reduce to 




*2 h ^2 ^2 “2 ^2 “2 


FJ . V (T' Oij dp. 


The equations of motion are linearized with respect to rotations about 
the y and z axes. That is, the angles and (i=l,2) are assumed 
small such that the following small angle approximations can be made: 


sin 0 . = 0 . 3 

ly ly 

cos 0 , = 1 , 

ly 


s in 0 , = 0 . , 

iz iz 

cos 0 . - 1 . 

xz 


In addition, the products ^of small angles and their rates ^ are neglected, 

i.e. , 07 , 0? , 0 . 0. , 0? , 0? , 0. 0. , 0. 0. , and 6. 0, are 

xy xz xy xz xy xz xy xz xy xz xy xz 

all assumed equal to zero. The coordinate transformation matrix between 

the inertial frame and body 1 (i = 1, 2) is developed from a 3-2-1 Euler 

angle sequence going from the inertial to the body i reference frame, 

that is , 


{x} = [-0. ] [-0. ] [-0. ] {x}. = T. {x}. 

I xz xy XX X X X 


where the nonlinear and linear.ized single axis rotation matrices 
defined as 



“ 


— 


~ 




1 

0 

0 


1 

0 

0 . 

ly 

“ 

0 

cos 6 , 

XX 

-sin0 . 

XX 

, [-0. ] = 
xy 

0 

1 

0 


0 

sin0 . 

XX 

COS0 . 
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-0. 

L ly 

0 

1 


are 


B-1 



NORTHRdP SERVICES, INC. 


TR-1400 


and 


[-6 




- 0 . 0 
xz 

1 0 . 
0 1 


Therefore 


T. = 
1 


" 1 

e, 

IZ 

-e. 

L ly 


e. 

ly 


sin6. - 0. 
XX xz 

cos6 . 

XX 

sin0 . 

XX 


COS0 . 
XX 


0 . COS0 , + 0 . sinO , 

xy XX xz XX 

-sin0 , 

XX 

COS0. 

XX 


From the kinematlcal relationship, 


T, = T. 0), , 

X XX 


the angular velocity of body i (i = 1, 2) is obtained as follows 




• 
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XX 
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XX 
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XX ^xz 

sin0 , 
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XX 

0. 
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XX 
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0. 
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xz 

COS0. - 0. 
XX xy 

sin0 . 

XX 


0. 
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0. 
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• The body fixed reference frames are assumed to be principal axes of 
inertia. Therefore, the body inertia matrix of body i (i = 1, 2) can 
be written as 


I. 

X 


"a. 0 o' 

X 

0 B. 0 . 

X 

0 0 c. 

L ij 


e Each body is assumed to be symmetrical about the body x-axis with 
respect to inertia, that is, B. = C, for i = 1, 2. Therefore, the 
inertia matrix of each body is^furtfrer reduced to 



0 

C. 

X 

0 
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• The center of mass of each body is assumed to lie on the x-axis . 

• The joint or hinge connecting the two bodies is also assumed to lie 
on the x-axis. Therefore, the joint vectors d^ and ^2 in their 
respective body reference frames can be expressed as 

, i = 1, 2. 



Based on the above simplifying assumptions, the matrix equation (equation 
(A-1)) given in Appendix A, describing the rotational motion of two coupled 
rigid bodies connected at a single point, expanded in scalar form can be 
written as 


A^ 6, = L, 


lx lx 


C*(0^ +6, 6. ) - K(0„ + 0„ 0- ) = L, cos 0, - L, sin 0, 

1^ ly Iz Ix^ 2y 2z 2x' ly lx Iz lx 


(A. - C. ) 0. 0, - y d, (d„ 0„ 0^ - d, 0, L ) 

1 1 Iz lx lx' 2x 2z 2x lx Iz lx 


m • 


C*(0. - 6. 0- ) - K(6„ “ 60 ) “ sin 0^ + cos 0^ 

1 Iz ly lx 2z 2y 2x^ ly lx Iz lx 


^2 ®2x ^2x 


• 9 


9 9 


+ (A, - C_) 0,,^ 0,^ + u d,^(d,,^ 0,„ 0,„ - d,„ 0,.. 0,. ) 


r ly lx 


lx'“2x '^2y 2x lx ly Ix^ 

(B-1) 


• • 


^2<®2y ®2z ^Zk> ‘ V ®lz " -^2y ®2x ’^2z ®2x 


- (A^ - Cj) + U ®iz ®ix^ 


# • 
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From 


then 




^2x ' hx + - ®lz> «lx - '®2y ‘ ®ly> ®lx hy 


(®2y - V> ®lx + <®2z - «lz> ®lx hz 


(B-2) 


Since 6, and 6, are assumed to be small angles, it is reasonable to assume 
ly Iz 

that the torques and will also be small in comparison to Conse- 

quently a first order approximation of equation (B-2) is 

^2x " ^Ix 

and the first and fourth equation of (B-1) becomes 


h ®lx " ^2 ®2x Hx' 


If in addition, it is assumed that 02x’ ’ then these two equa- 

tions decouple from the other four equations of (B-1) and can be solved separa- 
tely. Therefore, the values of ~ their derivatives can be 

regarded as known functions of time in the remaining four equations. Since 
the other four equations are linear in the other variables, they can be written 
in matrix form as 


J 0 + K 0 = L 


(B-3) 


where 
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K = 
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where 


L. = L. cos 6. - L. sin 0. 

ly xy XX xz xx 


L. = L. sin 6, + L. cos 6. 

xz xy XX xz XX 


i = 1, 2. 
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Appendix C 

SIMPLIFIED LINEARIZED EQUATIONS OF MOTION FOR TWO COUPLED 
RIGID BODIES IN TERMS OF NEW VARIABLES 


The two“-body equations Of motion as developed in Appendix B are in terms 

of the Euler angles ^2y’ ^Iz ’ ^2z* their rates and accelerations. 

Since the interest herein is concerned with the misalignment of the vehicle 

x-axes , it vjould be desirable to reexpres.s these two-body equations of motion 

in terms ofri=6„ - 0 t and 5 = 6_ - 0. , which provides a measure of 

2y ly 2z Iz 

misalignment directly. In addition, these two parameters are ones which can 
actually be sensed. 


In order to facilitate the above change of variables, it becomes necessary 
to define an additional set of variables, namely 


X = a 9 
Y = a 0 
Z -•= a 0 


lx 

ly 

Iz 


+ U-“>®2x 

+ (i-c)e2y 

+ a-“)02x- 


These variables are analogous to center of mass coordinates where a is a para- 
meter yet to be determined. Combining the above definitions for n, Y, and 
Z results in 
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and inverting gives 



The two-body equations of motion given by equation (B-3) can now be written as 
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? + (J b”^) ^ K b“'^ ^ 


-1 " 

= (J E ) L = F 


Expanding F results in 


[a C» + (1 ' c«)K] -[(I - a)CJ + aKl L 

[a C* + (1 - a)K] - [(1 - ci)c* + aK] 

(K - cj) + (K - C|) 

(K - C*) + (K - C|> £ 2 , 


In order to simplify F, a is chosen such that the first two elements of F 

A A A 

vanish whenever L, - L„ and L, - L„ . 

ly 2y Iz 2z 


Therefore, 


[a C* + (1 - a) K] = [(1 - a) CJ + aK] 


from which 


a = — : 

^1 


where C* - K and 


C* -K. 


Equation (C-2) can now be written as 
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MATHEMATICAL MODELS USED IN THE 2BODY SIMULATION 
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D.1 2BODY MODEL AND COORDINATE SYSTEMS 

The two coupled bodi:^s are modeled as two separate rigid bodies 
connected by a massless joint or hinge as illustrated in Figure D-1. The 
center of mass (CM) of each body is located relative to an inertial coordin- 
ate system by and ^ 2 * Although the simulation assumes that initially 
(at time = 0), = Rj^ = 0, there is no loss of generality. The position 

vectors and ^2 locate the joint or hinge point of each body relative to 
the center of mass. 



Figure D-1. TWO COUPLED RIGID BODIES INERTIALLY REFERENCED 

The relationship between the inertial and body coordinate frames is 
as follows: 

{X}j = {X>3^ 

{X}^ = T 2 {X>2 

(X)^ - T2{X>2 = Tj {X)2 = T'CX)^ 

The transformation T^, (i=l,2>, is devel.oped from 3-2-1 Euler angle sequence 
going from the inertial to the body i reference frame. 

D-1 


ft 





D.2 JOINT MODEL/CONSTRAINT AND ALIGNMENT TORQUE 

The joint connecting the two rigid bodies is treated as an ideal docking 
mechanism during the soft dock regime (the bodies are physically attached, 
but not yet hard docked (rigidly attached), angular motion can exist between 
the two vehicles, and the target vehicle could be spinning about an axis 
close to or coinciding with the docking port axis) . 


The connection or joint is modeled as a massless spring and damper. 
From Figure D-2, it is seen that 



Figure D-2. TWO COUPLED RIGID BODIES WITH IDEAL JOINT 


A “ ^2 ” — 1 — 2 ~ — 1 

and 

A “ ^2 “ ^2 ^ ^2 “ M.]l ^ 

where it is assumed that << |Aj^|> i ~ 1,2. 
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Initially (at time = 0), however, 
d = d = 0 

since the constraint relations 


and 


^2 * + 4 ■ 4 

R„ = R-, + (D, X d^ - w- X d 


are used. The joint structure provides a physical constraint force and 
torque given by 


F_ = K d + C d 


and 


L = d X F 


where K and C represent the structural flexibility and damping coefficients. 
To damp the scissoring-type angular motion between the vehicles, an alignment 
torque is provided, i.e.. 


V + C 


A 


V 


where and are spring and damper coefficients and ^ is an error vector 
based on the amount of angular misalignment between the docking port axes. 
The program assumes the x-axes of the vehicles are to be aligned. Thus 


V = ii X i2 


where ij is the x-axis unit vector of body j, j = 1,2. Since 
then V referenced to the body 1 coordinate frame is 


V = i^ T' ., = 
— 1 jl 


0 

-T' 

31 

T' 

21 


where the cross product matrix operator of a general vector 
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is formed by 



where the elements of v can easily be determined from 
T' = T' 0)2 - 0)^ T’ . 

To determine values for K and C., the alignment torque is assumed of 

A A 

the form (also a typical rate-position feedback control law) 

L = I 6 = - 0 - e 

where 

I - an effective inertia of the sytem 
9 - a relative angular displacement 
9 - a relative angular rate. 

Rearranging terms results in 

.. C K 

O+~9+^0 = O 

which is recognized as the equation of motion for a linear harmonic oscillator 
where 
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2 

0 ) 


Thus 

Ka = I I 

and 

CA=2CIto = 4TT I?f 

where 

f - undamped natural frequency 
C - non-dimensional damping factor. 


The undamped natural frequency associated with alignment should usually 
be chosen to be lower than the precessional frequency of the spinning body. 
The effective inertia can be approximated by that of the reduced inertia of 
the system, i.e.. 


I = 


h "2 
h + ^2 


where I, is the inertia of body j, j = 1, 2. 

1 

D.3 DESPIN TORQUER MODEL 

When one body is spinning, a torquer is available for despinning. The 
despin torque can be applied from either body (reflecting actual hardware 
considerations), but such application is on only the x-axis . For high 
relative spin rates, the despin torque is limi.ted, providing a constant 
output, whereas for low relative spin rates the torquer output is linear as 
seen by the following equations and illustrated in Figure D-3. 


Ljjj = sgnCtOj^ - Tlim , -A» > 

Sx " ^ “2x - “lx " 


The despin torque may be either excluded from the simulation, turned on 
at time = 0, or activated on an energy condition. The energy condition might 


D-5 



NORTHROP SERVICES, INC. 


TR-1400 



^DX 

- 

/ 

- Aco 

/ Acj 

t , 


7 


/. 

1 

k —T 


'lim 


Figure D-3. DESPIN TORQUER OUTPUT MODEL 


be desirable if large initial misalignment conditions were used. Only after 
a highly transient condition damped out would it be appropriate to despin. 


D.4 BODY 1 ATTITUDE CONTROL 

If attitude control of body 1 (usually throught of as the chase vehicle, 
whereas body 2 is usually thought of as the target and possibly spinning vehicle) 
is desired, an ideal control system is available. The angular rates of body 1 
are continuously set to zero to maintain a fixed vehicle attitude; hox<rever, the 
vehicle is" still free to translate. 
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